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TITLE 

NUMBER OF WORDS 

TIME 

TYPE 
PURPOSE 

TEMPqitARY STORAGE 
PARAMETERS 



D^CRth?IO!r 



Matrix Multiplication (Closed) (DOI or SADOl) 

77 words 

[(i x k) (2d A + 1) + M(k + j) (2d B + 1) + (i x j) 2 k] ms 

(See note 1.) 

Closed 

To multiply two matrices A»B » C, retaining the product 

matrix C in the storage. 

0, 1, I) - k, b - k + 1, . * -, t . 1 

S3? 00 JF 00 bF 

j » number of columns in B 
b = location of the first element of B 
The routine is entered at word from one of two entries; 



a) 



p 


kO r 
50 p 


(2) 


P 


NO r 
50 p 




2k to routine 




2k to routin 



r « location of the first element of C 
Entry (l) assumes neither A ■ ( a *v) n or B ■ (\*) is 
contained in the storage. Entry (2) assumes B is contained in the storage but 
A is noi* In either case, before proceeding vith the multiplication the entire 
matrix gj must be contained in the storage, then matrix A is read in one row at 
a time**tne (k + l)st row replacing the k th row— forming C. The elements of a 
row of A immediately precede B, i.e*, a u « H(b - l), a » w(b * 2), 

* * * a ik - *^ * ^ ' Heme ^e "total number of location? required for matrix 
element* is k + (k x j ) + (i x j ) , 

, Two checks are incorporated in this routine. One check teste 
to see that -1 < c., < 1. Whenever overflow occurs (i»e. 1 < c < -1) the com- 
puter will stop on an FF order in word 4lL. The other Check tests to see that 
all rows of a matrix being read into the storage contain the same number of elements , 
j^ any ^q TCfW8 of fl j^^pj^ con tain a different numbe^f of elements the computer 
will stop on a FF order in word 2%. 



-2 - 



METHOD OF USE To avoid overflow, & ± ^ and b must be scaled so that 
|c | = £ la^ b | < 1. This will be true if ja^l < x/n and [b^ |< x/n 



2 . 
for x < n. 

Prepare Matrix Tape A by punching the elements (from left 
to right) of the first row of A followed by the second row, etc., as follows: 

a 00* a 01' — a 0k' N ' a 10' a ll* — a lk' N * ' a iG* a il' — *ik' *' J ' ^ eh 

element must be preceded by a plus or minus (K or S) sign which can be followed 
by up to 11 decimal digits* the last element of each row must be followed by Nf 
and the last element of the matrix must be followed by J in addition to If. 

If B is to be read from tape, Ifetrix Tape B should be pre- 
pared in a similar manner. If B is contained in the storage, i.e., Entry (2) 
is used, the elements of B must be stored in the same sequence in which they 
would appear on tape, but the row and matrix ending symbols H and J should not 

be used. 

When Entry (l) is used the computer will stop twice, first 
on a 2k order and then on a 20 order. At the first stop Matrix Tape B should be 
placed in the reader and at the second stop Matrix Tape A should be placed in 
the reader. After all stops computer operation should be resumed by moving the 
black switch. 

H0TB3 (1) A is an (i x k) matrix 

B is a (k x j) matrix 

d. is the number of decimal digits in an element 
A 

of matrix and likewise for d_. 

M ■ 1 if Entry (l) is used. 
* if Entry (2) is used. 
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lgr 



LOCATION 



3 
4 

5 
■6 

7 
8 

9 
10 
11 
12 

13 
14 

15 
16 

17 



ORDER 



00 K(ML1) 

K5 F 

42 76L 
10 20F 
42 4% 
36 % 
F5 25L 
42 25L 
L5 671^ 
4o i8l 
L5 65L 
40 30L 
L5 63L 
40 21L 
L5 33L 
42 13L 
00 2 OF 
46 13L 
40 F 
81 4F 
LO 6lL 
32 15L 
L4 6lL 
JO 58L 
10 4F 
75 F 
00 4F 
L4 F 
40 F - 
50 F 
7J 60L 
22 8L 
LO 58L 
32 22L 
L5 60L 
40 F 
L5 13L 
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Mil 



Set address of c 

Test for N or 4 



U 



_ Prepare to read Tape A 



- Set address of a 



11 



. Form a., or "b, . 
^ ik kj 



Test for K(+) 



Q 



LOCATION 


ORDER 
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1 l8 




LU 59L j 






19 


. 


J+2 19LJ 
*+0 13L 




— Set address of a,, or b, . f 

1 




1+1 F 




* 

! 


20 


F5 62L 
kO 62L 




1 

Count every row | 


21 
22 




^5 2l+Ll 
1+0 21+Lj 

26 9L 




1 
Count first row 




LO 58L 






23 


36 25L 
LI 60L 




Test for S(-) 

- ■ 1 


2k 


26 17L 




i 

? 




00 P 




1 
,1 counter j 


25 


LO 58L 




i 




32 68L 




Test for N 


26 


L5 22L 








kO 21L 






27 


L5 62L 
LO 2*+L 




. 


28 


1+0 F 
L3 F 




— Test number of elements in each row 


29 


36 30L 




! 

1 




FF F 




| 


30 


Po 62L] 




1 


31 


IM- 33L 
00 20F 






32 


k6 68L " 
1+1 IF 




1 




JO IF 
S5 F 




! 

ij 


3^ 


50 F 
lh F 




4 

j 


35 


kO F 








32 38L 







LOCATION 


ORDER 
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36 


L5 IF 
32 1+1L 








37 


L*+ f 
36 1+1L 








38 


26 1+2L 
L5 IF 








39 


36 1+OL 
22 1+1L 






- Form c, , and test for overflow 


1+0 


L4 f 

36 1*24 








in 


FF F 
LU f 








1+2 


1+0 IF 
L5 57L 








h3 


L*+ 3^L 
1+0 5^L 


* 






kk 


LO 68L 
32 33L 








^ 


L5 IF 
1+0 F 








k6 


F5 1+5L 
1+0 1+5L 






hi 


F5 62L 
1+0 62L 






_. Has row i of C been completed? 


1+8 


LO 56L 
32 51L 








V9 


L5 55L 
1+6 3^L 






50 


Ll+ 62L - 








1+2 3^L 


. 




51 


22 32L 




Return to calculate c. , , 




V? 55L 


. 




52 


$+0 31+L 






. 


i+3 62L 






53 ! 


£2 6L 




Calculate c. , , 




00 59L 




b , 



LOCATION 



ORDER 



FOTES 



PAGE k 



M II 



5^ 


iAf 




l<-0 F 


55 


50 F 




7*+F 


56 


26 17L 




00 F 


57 


LL 4095F 




00 F 


58 


00 F 




00 IF 


59 


L5 3F 




1+2 3^L 


60 


1+2 55L 




1+2 33L 


61 


LO 58L 




1+2 13L 


62 


1+2 5^L 




10 20F 


65 


1+2 56L 




1+2 57L 


61+ 


L5 13L 




00 2 OF 


65 


1+6 5^L 




1+6 13 L 


66 


1+6 3^L 




h6 55L 


67 


22 67L 




L5 531/ 


6a 


22 1011+F 




26 59L 




26 IN 


59 


00 IF 




00 IF 


60 


ON 3276F 




NN 3277F 


61 


00 F 




00 10F 



- Constants 



- Interlude 
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LOCATION 



52 



63 
64 

65 
66 
61 
68 
69 
TO 

Tl 
72 

73 

74 

7=i 



ORDER 

26 17L 
00 F 
F5 24L 
40 24L 
43 62L 
26 9L 
43 62L 
LI 2^4-L 
L4 59L 

42 19L 
LO 59L 
1+2 19L 
50 F 

43 24L 
20 2L 
L5 5^L 
40 13L 
L5 66L 
Ho 18L 
L5 63L 
40 21L 
1+3 24L 
L5 22L 
L4 6lL 
42 25L 
L5 64L. 
40 30L 
L5 55L 
40 3^L 
22 F 



NOTES 



PAGE 5 



M-L-" 



Constants 



Reset initial conditions and address ss 



Link 



